Year Two of Archaeofaunal Analysis at the Haynie Site (5MT1905)

report
zooarchaeology
taphonomy
machine learning
Authors
Affiliations

Jonathan Dombrosky

Crow Canyon Archaeological Center

Department of Anthropology, University of New Mexico

Ahna Feldstein

Department of Archaeology, Boston University

Published

February 8, 2024

1 Introduction

The Crow Canyon Archaeological Center initiated the Northern Chaco Outliers Project (NCOP) in 2016 (Ryan 2016), and fieldwork began in earnest in 2017 (Simon et al. 2017). The NCOP’s main purpose is to study the Lakeview Community, which is a collection of four neighboring Chaco period Great Houses located in southwestern Colorado. The four Great Houses are subsumed by three archaeological sites: the Haynie site (5MT1905), Ida Jean (5MT4126), and Wallace Ruin (5MT6970).

Excavations at the Haynie site (5MT1905) have comprised the vast majority of NCOP-related fieldwork over the last six years. Haynie is a complex, multicomponent ancestral Pueblo archaeological site intermittently occupied between roughly A.D. 800 and 1175 (Throgmorton et al. 2022, 2023). It includes two Great Houses, a large Pueblo I-II roomblock, and numerous other structures and contexts dating to the Pueblo I and II time periods.

Examining archaeofaunal remains from the Haynie site is one important part of the Northern Chaco Outliers Project. In the first year of archaeofuanal analysis, Dombrosky and Gilmore (2023) reported three different findings worth summarizing here. First, the Haynie fauna has a low rate of identifiability associated with a high amount of fragmentation. Preliminary predictive models indicated that thick cortical bone (from medium mammals and larger) fragments are highly associated with unidentifiabiltiy. However, those preliminary models did not take into account whether unidentifiable remains are associated with disturbed contexts, which is a conspicuous feature of the Haynie site. Second, identification rate estimates suggested that new identification types were still quickly added as we progressed through the analysis. Those results indicated that our estimations of taxonomic diversity were still incomplete, but that the Haynie site was still likely not notable in terms of taxonomic diversity compared to other sites in the central Mesa Verde region. Finally, the first year of archaeofaunal analysis also focused on two different rare taxa recovered from the Haynie site: wolf (Canis lupus) and bison (Bison bison). We developed a highly accurate model used to classify wolf, coyote (Canis latrans), fox (Urocyon spp. and Vulpes spp.), and domesticated dog (Canis familiaris) from mandibular measurements using a massive reference database provided by Welker et al. (2021). The model we developed was 100% confident in identifying the mandible specimen from Haynie as wolf. Additionally, we illustrated how we identified bison remains and discussed the clear importance of medium-to-large game procurement at the site.

In this report, we revisit and expand some of these previous topics all while reporting on a year’s worth of new identifications. In the identification section (Section 3), we repeat some of the information presented before, but have provided a new subsection that focuses on summarizing the contexts that dogs (genus Canis) have been recovered from in the central Mesa Verde region and how the Haynie site compares. The new section called Taxonomic Abundance by Study Unit (Section 4.1) works to identify excavation contexts that are most likely disturbed, and the taphonomy section (Section 4) incorporates this new information into a new model to tease apart factors driving unidentifiablity. Finally, we re-estimate identification rates at the site to gauge taxonomic representativeness and compare these estimates to all of the faunal assemblages in the Crow Canyon database. Our larger goal is to document the quality of zooarchaeological data produced from May 2023 to August 2023; we firmly believe that a focus on data quality is essential to valid archaeological interpretation. Data quality underlies almost every facet of archaeological research.

2 Materials and Methods

There are three analysts associated with the data described here: Jonathan Dombrosky, Eric Gilmore, and Ahna Feldstein. Jonathan Dombrosky has approximately 12 years of experience with archaeofaunal analysis, and he analyzed new specimens associated with this report from May 2022 to August 2023. Eric Gilmore had approximately 3 years of experience with archaeofaunal analysis, and he analyzed specimens as a Crow Canyon Zooarchaeology Intern in the summer of 2022. Ahna Feldstein has 4 years of experience with faunal analysis and was a Crow Canyon Zooarchaeology Intern from May 2023 to July 2023. Jonathan, Eric, and Ahna are respectively referred to as Analyst 1, Analyst 2, and Analyst 3 in all subsequent interanalyst comparisons.

The main comparative collection used for the analysis of the Haynie site archaeofauna is housed in Crow Canyon Archaeological Center’s Laboratory. Three cottontail (Sylvilagus spp.) specimens and one black-tailed jackrabbit (Lepus californicus) specimen were loaned from the Laboratory of Zooarchaeology at the University of North Texas. We took specimens that were difficult-to-identify to the Museum of Southwestern Biology’s Division of Mammals and Division of Birds located in Albuquerque, New Mexico. Various osteological guides, manuals, atlases, and keys aided identification. Several publications assisted with the identification of mammal remains (Adams and Crabtree 2012; Chavez 2008; Gilbert 1980; Hillson 1986, 1996; Jacobson 2003; Olsen 1964; Smart 2009). Bird remains were identified with other works (Cohen and Serjeantson 1996; Gilbert et al. 1981; Hargrave and Emslie 1979; Olsen 1979), as well as fish, amphibian, and reptile bones (Olsen 1964). Some consulted references helped identify both avian and mammalian remains (Broughton and Miller 2016; Elbroch 2006). Yet other works verified nonhuman from human remains (Baker et al. 2005; France 2009; White et al. 2012).

We followed identification protocols explicitly designed to enhance data quality (Driver 1992, 2011; Wolverton 2013; Wolverton and Nagaoka 2018), and used the coding system by Driver (2006). Briefly, analysts adopted a conservative approach to identifying zooarchaeological specimens at the Haynie site. It is an almost impossible task for analysts to understand how all diagnostic skeletal criteria change through time, among species, within different age classes, between sex, and across geographic areas on a fragment-by-fragment basis. It has been argued that identifications become less taxonomically specific when analysts have more experience, greater access to diverse comparative materials, and a specific focus on data quality (Gobalet 2001; Lyman 2002, 2019; Wolverton and Nagaoka 2018). This lack of taxonomic specificity likely increases identification accuracy in situations where assemblages contain an abundance of fragmented remains from closely related taxa.

We use the Number of Identified Specimens (NISP) to report taxonomic abundance, and this quantitative unit is a tally of all archaeofaunal specimens within a given taxonomic classification. NISP is the most basic quantitative unit from which most others are derived, such as the Minimum Number of Individuals (MNI). NISP is preferred because it is often highly correlated with measures like MNI. It is also devoid of errors in additive calculation that plague minimum number units (Grayson 1984; Lyman 2008). We also rely on a non-standard unit called Unique Identification Types (UITs) to estimate taxonomic diversity (Section 5). In the last report, we called this unit the Number of Unique Identifications (NUIDs). The way this unit is calculated is exactly the same. We have simply changed the name after conversations with colleagues.

All statistical analyses and figures were produced with R version 4.3.1 (R Core Team 2023). Our statistical analyses are structured with tidyverse packages and syntax (Wickham et al. 2019). All graphs were produced with ggplot2 (Wickham 2010). We built predictive taxonomic and taphonomic models using a supervised learning workflow (Hastie et al. 2009; James et al. 2013; Kuhn and Johnson 2013); this included using the tidymodels metapackage to split our data and implement basic model features (Kuhn and Silge 2022). We rely on logistic regression and random forest as the engines for our predictive models. Logistic regression is a powerful modeling engine designed for binary classification (Kuhn and Johnson 2013, 282). Random forest models predict classes based on one or more variables through a large assortment (a forest) of small decision trees (Kuhn and Johnson 2013, 198). Random forest models sometimes outperform more traditional predictive models (Cole et al. 2022).

3 Identified Taxa

There are 2191 identifiable and 5640 unidentifiable specimens so far in the Haynie archaeofaunal assemblage. We have added 590 identifiable and 1359 unidentifiable specimens—for a total of 1949 specimens—since our last report. Our identification rate is still low at 27.98%, which has practically gone unchanged. Of all of the large faunal assemblages that Crow Canyon has analyzed over the years, Haynie currently has the lowest identification rate (Figure 1). Analyst 1 analyzed 62.47% of the assemblage, Analyst 2 analyzed 20.6%, and Analyst 3 analyzed 16.93%. There are 4 classes of animals present: Mammalia (mammals), Aves (birds), Actinopterygii (ray-finned fishes), and Reptilia (reptiles). There are 4 orders of mammals present, 7 orders of birds, 1 order of ray-finned fishes, and 1 order of reptile. In total, we used 61 identification types (Figure 2).

Figure 1: The Haynie site (highlighted in light blue) currently has the lowest identifiablity rate of Crow Canyon faunal assemblages that exceed 1000 Number of Identified Specimens (NISP). Percent identifiable is the ratio of identifiable specimens to all specimens in a faunal assemblage. Number of Specimen (NSP) values for each site are located within each bar. It currently ranks as the 48th lowest identifability rate in the enitre Crow Canyon database among 52 sites.

Figure 2: Relative taxonomic abundance at the Haynie site using percent Number of Identified Specimens (%NISP). Raw NISP values are reported next to each bar. Total NISP = 2191.

3.1 Mammalia (n = 1964)

3.1.1 Lagomorpha (n = 761)

Cottontails (Sylvilagus spp.) and jackrabbits (Lepus spp.) are the main taxa in the order Lagomorpha at the Haynie site. Lagomorphs are currently the most abundant animals identified, comprising 34.73% of the identified specimens. It has been hypothesized that populations of larger-bodied jackrabbits decreased through time in the central Mesa Verde region, and that this decrease is likely due to human overhunting (Driver 2002; Ellyson 2014). Thus, the ratio of cottontails to jackrabbits—commonly referred to as the Lagomorph Index—is a basic quantitative unit of general interest in the area, and in western North America in general (Driver and Woiderski 2008). The Lagomorph Index at the Haynie site is 0.68, which indicates a fairly equal relationship between the abundance of cottontails and jackrabbits. Haynie currently has the lowest Lagomorph Index of all the large (greater than 1000 NISP), central Mesa Verde faunal assemblages in the Crow Canyon database (Figure 3). This number is similar to the Lagomorph Index from some Pueblo II components of Shields Pueblo (5MT3807). Shields serves as an important point of comparison—here and in subsequent analyses—considering that it and Haynie both have similar site features (i.e., Great Houses) and general occupation histories (Rawlings 2006).

Figure 3: The Lagomorph Index across all faunal assemblages exceeding 1000 Number of Identified Specimens (NISP) in the Crow Canyon database. Currently, the Haynie site has the lowest Lagomorph Index of all the central Mesa Verde assemblages (Tsama Pueblo is a Pueblo IV Northern Rio Grande site).

Garden hunting is an important subsistence practice to consider at the Haynie site, and the moderate Lagomorph Index value is also interesting in this regard. Researchers argue that higher ratios of cottontails to jackrabbits indicates a higher reliance on garden hunting in the central Mesa Verde region (Driver 2011). What could the more even relationship of cottontail and jackrabbit abundance indicate about garden hunting at the Haynie site? Future work using stable isotope analysis (sensu Dombrosky et al. 2023) might help shed light on whether garden hunting was prevalent at the site, and it also might help describe the relationship between the Lagomorph Index and garden hunting in general.

3.1.2 Small Mammal (n = 339)

This identification group includes those mammals jackrabbit size and smaller. This non-standard identification includes all small mammal specimens lacking morphological features required for more specific taxonomic levels. Consequentially, it has a high likelihood of incorporating many lagomorph specimens since they are exceedingly abundant in southwestern archaeofaunas. This identification group will be incorporated into final comparisons of size-based abundance indices through time, between areas of the Haynie site, or between sites. The goal will be to gauge the impact this identification group has on final interpretations of the Lagomorph Index.

3.1.3 Artiodactyla (n = 340)

Even-toed hoofed animals make up the order Artiodactyla, and they are currently 15.52% of identified specimens at Haynie. Substantial zooarchaeological evidence indicates artiodactyls—mostly mule deer (Odocoileus hemionus), pronghorn (Antilocapra americana), and bighorn sheep (Ovis canadensis)—were severely overhunted in the central Mesa Verde region (Badenhorst and Driver 2009). Some models suggest that deer populations might have been so low that any observable deer was immediately hunted after A.D. 1000 (Bocinsky et al. 2012). For this reason, the Artiodactyl Index is another basic quantitative unit of general interest in the region and in the larger northern U.S. Southwest. It measures the ratio of large-bodied artiodactyls to small-bodied lagomorphs in an archaeofaunal assemblage (Broughton et al. 2011). The current Artiodactyl Index at Haynie is 0.31. This value may seem somewhat low, but in reality it is quite high. The Haynie site has the highest Artiodactyl Index of all the large (exceeding 1000 NISP) faunal assemblages in the Crow Canyon Research Database, the only exception is a Pueblo IV site from the Northern Rio Grande called Tsama Pueblo (Figure 4).

Figure 4: The Artiodactyl Index across all faunal assemblages exceeding 1000 Number of Identified Specimens (NISP) in the Crow Canyon database. Currently, the Haynie site has the highest Artiodactyl Index of all the central Mesa Verde assemblages (Tsama Pueblo is a Pueblo IV Northern Rio Grande site).

Large artiodactyls—elk (Cervus canadensis) and bison (Bison bison)—are also notable parts of the Haynie archaeofaunal assemblage. These specimens comprise 6.47% of artiodactyls. There are currently 4 bison specimens, and Dombrosky and Gilmore (2023) reported on these specimens in detail. Importantly, these specimens do not overlap with each other and represent an MNI of 1.

We radiocarbon dated five bison specimens to help determine whether more than one bison individual is present at the Haynie site (Figure 5). Four specimens were found in somewhat similar upper layers of a trench (PDs 480, 489, and 505). One specimens (PD 1265) was found in lower layers of this trench within a feature. The conventional radiocarbon age and error of the specimen from the lower feature is identical to two specimens from the upper layers, while two separate specimens from these upper levels have somewhat different calibrated ages. Altogether, there is almost complete overlap in the 95% calibrated date ranges. Radiocarbon dating has not confirmed whether there are separate individuals or not, but what is apparent is that either deposition occurred rapidly in this general area of the site or that stratigraphic mixing is more pernicious than originally anticipated.

Figure 5: Five calibrated radiocarbon dates from bison specimens at the Haynie site. All specimens calibrated using IntCal20. C:N atomic values of collagen were all within 2.9 and 3.6, which is indicative of well-preserved collagen (Ambrose 1990).

3.1.4 Medium Mammal (n = 202)

Mammals larger than a jackrabbit and up to deer size are considered medium mammals. There is a high probability this identification group incorporates many artiodactyl specimens since they are one of the most common medium mammals in southwestern archaeofaunas. It should be incorporated in sensitivity analyses that rigorously assess final conclusions relying on size-based abundance index values, such as the Artiodactyl Index.

3.1.5 Rodentia (n = 187)

Rodents are 8.53% of identified specimens, which is a small component of the overall assemblage. The ratio of sciurids to rodents can help gauge the impact of intrusive species. Members of the squirrel family (Sciuridae) are notorious intruders of archaeological deposits, especially prairie dogs (Cynomys spp.). The majority of prairie dog skeletal fragments are identified to the family level, as it is extremely difficult to skeletally distinguish prairie dogs from ground squirrels. This means that the Sciuridae identification has the highest potential for accumulating prairie dog specimens. And, indeed, we can confirm that the vast majority of sciurid specimens compare favorably to prairie dogs. Sciurids make up 39.04% (n = 73) of the rodent assemblage. Given that rodents comprise a small portion of the overall assemblage, it does not appear that rodent intrusions pose a significant problem for interpretation at the site. We will, however, closely track the taxonomic composition of rodents as the Northern Chaco Outliers project progresses. Importantly, rodents could have been actively hunted (Badenhorst et al. 2023). One way to disentangle intrusive from non-intrusive rodents is through radiocarbon dating (Guiry et al. 2021). This could prove useful if prairie dog specimens preclude clear interpretation of the Haynie archaeofaunal assemblage in the future.

Large rodents, those larger than a woodrat (Neotoma spp.), comprise 18.72% of the total rodents currently identified. Beaver (Castor canadensis) and porcupine (Erethizon dorsatum) specimens are the most notable. Beavers prefer aquatic habitats while porcupines take refuge in trees along active floodplains (Baker and Hill 2003; Roze and Ilse 2003). The presence of these species might suggest hunting activities that were focused in riparian habitat close to the site.

3.1.6 Carnivora (n = 119)

Carnivores are a small portion of currently identified specimens at the Haynie site (5.43%), but a number of canid specimens were discussed in detail in last year’s report (Dombrosky and Gilmore 2023). There are currently 76 specimens from the family Canidae in the Haynie site fauna. Many of the specimens come from an articulated domestic dog offering and one specimen is a wolf mandible. We provide a contextual analysis of Canis spp. deposition across all Crow Canyon projects considering the importance of this group of animals to the Northern Chaco Outliers Project.

3.1.6.1 Spatial Context of Canis spp. Specimens

Domestic dogs and their wild counterparts filled many roles within the Ancestral Pueblo world; dogs were used for companionship, raw materials, subsistence, and in ritual practices such as dedicatory offerings in burials, and artistic depictions (Monagle and Jones 2020). Recent investigations from Arroyo Hondo Pueblo suggest that, contrary to our modern interpretations of domestication, no such dichotomy existed between domestic and wild canids in Pueblo society (Monagle et al. 2018). Coyotes and dogs often swapped or shared roles within their communities, and both species were included in instances of ritual deposition (Monagle et al. 2018).

It is clear that canid remains offer great insight into past human lifeways (Semanko and Ramos 2022), and they hold an especially important space within the cultural landscape of the central Mesa Verde region. Aside from the turkey (Meleagris spp.), dogs are commonly cited as the only domesticated animal within this region prior to Spanish colonization. Ethnographic research on historic and modern Pueblo communities suggests that canids fulfilled specific roles as protectors, warriors, companions, hunting assistants, spirits, messengers, pests, trash-eaters, food, witches, and sometimes even substitutes for humans (Monagle and Jones 2020). The time-depth and spatial patterning of canine deposition is highly diverse in this region (Hill 2000). The Haynie faunal assemblage contains specimens attributed to domestic dogs, coyotes (Canis latrans), and wolves (Canis lupus), meaning that human-canid relationships are especially complex at this site.

For estimating the relative abundance of canid specimens at Haynie, NISP was chosen as our basic unit of analysis. As an additive measurement of taxonomic abundance, NISP avoids those issues of aggregation which often accompany derived units such as MNI (Minimum Number of Individuals). NISP and MNI nonetheless often share a strong correlation which may be examined to address issues of fragmentation and interdependence (Grayson and Frey 2004).

There are currently 70 specimens from the genus Canis in the Haynie archaeofaunal assemblage. The genus Canis includes domestic dogs, coyotes, and wolves. Taxonomic classification beyond the genus level is notoriously difficult based on post-cranial morphology, but select specimens from Haynie have been identified to the species level. Taxonomic identifications based on unique Provenience Designation (PD) are represented by Figure 7.

Canid specimens are relatively rare across the entire Crow Canyon archaeofaunal database, contributing less than 8% of total NISP at all sites (Figure 6). Canid remains at the Haynie site comprise a mere 0.89% to the total raw NISP of 7831. Given the established role of dogs, coyotes, and wolves within Pueblo society, the overall scarcity of canine remains indicates a selective process of canine deposition.

Figure 6: Relative abundance of all specimens identified to the genus Canis across the Crow Canyon faunal database. Raw total NISP values per site are reported next to each bar.

An initial spatial analysis of canine deposition at Haynie helps untangle the complex relationship between humans and dogs within Ancestral Pueblo society. By observing the abundance and spatial distribution of their remains, we begin to understand key associations between canine individuals and the built environment. A canid’s transition from life to death is preserved in its deposition, allowing us to glean information on how these individuals were felt, seen, and treated by their human counterparts. Continued contextual analysis of canid remains offers greater potential to understand how the Haynie site may adhere to or diverge from local trends of animal management within the broader Mesa Verde region.

To evaluate the spatial distribution of canid remains, we first pursue analysis based on provenience designations (PD) for each specimen at Haynie. According to the Crow Canyon Field Manual (Center (2001)), an artifact’s PD is assigned based on its mode of deposition and stratigraphic position.

The Haynie assemblage contains two important instances of canid deposition from the first year of archaeofaunal analysis (Dombrosky and Gilmore 2023). The first specimen, a wolf mandible, was recovered from the fill between two floors of a surface room (Throgmorton et al. 2022). The mandible from year one is currently the only specimen identified to Canis lupus (see Figure 7), confirmed through the application of predictive modeling with mandible biometrics.

The second set of specimens belong to a single individual, likely a small domestic dog. 56 of the total 70 canine specimens at Haynie belong to this domestic dog, all sharing one distinct PD. The dog’s upper half was uncovered from the floor of a pitstructure, and its lower half was deposited outside of the unit (Throgmorton et al. 2022).

During this year’s archaeofaunal analysis at Haynie, another domestic dog mandible was added to the database. The mandible was identified to Canis familiaris based upon mandibular body width, small overall size, and lack of elongation compared to coyote, wolf, and fox (Vulpes vulpes).

Specimens classified as Canis spp. are designated to the genus Canis, but cannot be identified to domestic dog, coyote, or wolf based on morphological characteristics. Specimens in this category may originate from any of the aforementioned Canis species. Coyote remains are likely represented by the Canis spp. category, though no specimens from Haynie have currently been identified to Canis latrans.

Figure 7: Number of Identified Specimens (NISP) identified to the genus Canis from the Haynie Site (5MT1905). The genus Canis includes domestic dogs (Canis familiaris), coyotes (Canis latrans), and wolves (Canis lupus). A majority of specimens at Haynie are attributed to one instance of domestic dog deposition, recovered from the floor of a pitstructure (Throgmorton et al. 2022; Dombrosky and Gilmore 2023). The Haynie faunal asseblage also features a rare instance of wolf deposition, and a domestic dog mandible added during the most recent year of faunal analysis.

To begin to examine depositional trends at Haynie, we employ descriptive categories established by Erica Hill to understand general trends across all sites within the Crow Canyon faunal database. The three main descriptive categories of canid deposition include ceremonial trash, dedicatory offering, and simple interment (Hill 2000). For this analysis, we apply the term deposition as opposed to interment to avoid assumptions of intentionality. We also change ceremonial trash to ceremonial discard in the following analysis, in adherence with historic and current Pueblo perspectives. We suggest the future use of ceremonial discard as opposed to Hill’s original terminology, the problematic nature of which is discussed below.

Dedicatory offerings dominate, with the highest NISP counts coming from surface burials (Figure 8). However, specimen interdependence is driving this trend: the articulated remains of a single dog comprise the majority of identified canid specimens at Haynie (Throgmorton et al. 2022; Dombrosky and Gilmore 2023). The small dog’s remains contribute to an inflated NISP count and MNI of 1. Greater overall NISP counts within surface burials may not, therefore, suggest the presence of more dogs; but instead the careful deposition and preservation of select individuals.

The domestic dog’s skeletal articulation and placement adhere to the key characteristics of dedicatory offerings, as outlined by Hill (2000). Dedicatory offerings feature deposition within ritual spaces, often representing the ceremonial closure of architectural features such as kivas. Similar deposition practices have been observed at sites within Mancos Canyon (Monagle and Jones 2020), and from Basketmaker III/Pueblo I period sites in Dolores, Colorado (Hill 2000), illuminating a broader pattern of ritual deposition across the Mesa Verde region.

The deposition of a singular domestic dog mandible in midden reflects our expectations for ceremonial discard. Ceremonial discard involves the ritual dispatch of an individual and their subsequent burial outside habitation areas, usually in a sacred space or shrine. Hill (2000) states that dogs suspected of witchcraft may appear as ceremonial trash. Midden deposits are often regarded as sacred spaces in Pueblo communities (Fladd et al. 2021), and the deposition of objects within middens is considered an important rite-of-passage between the lived and spiritual worlds (Walker and Berryman 2023). This mandible is the most recently-identified domestic dog specimen, and the lone example of ceremonial discard in the Haynie assemblage.

Figure 8: Spatial distribution of canid remains at the Haynie site, according to skeletal part frequency and deposition type.

Hill’s descriptive categories offer an avenue for distinguishing patterns within a complex system of human-animal interaction. To gauge how Haynie compares to other sites within the Mesa Verde region, the relative abundance of canid remains per depositional category is assessed across the entire Crow Canyon database (Figure 9).

Across all sites in the database dedicatory offerings are most abundant, followed closely by simple deposition. Simple deposition denotes instances of canid deposition lacking the characteristics of both ceremonial discard and dedicatory offerings. Remains within this category may indicate expedient disposal, focused on the discard of sick or burdensome canines (Hill 2000). Simple deposition may also contain particular forms of deposition which involve ritual behavior, but do not adhere to Hill’s classifications (Hill 2000).

NISP values for dedicatory offerings and simple deposition do not differ substantially, though instances of ceremonial discard remain scarce. Less than 2% of canid specimens across all sites fall into the category of ceremonial discard (Figure 9). Limited occurrences of ceremonial discard may imply the limited deposition of canine remains in middens, or their use in depositions outside the scope of Hill’s categories. The wolf mandible deposited between two floors is included in simple deposition, and is further discussed below.

Despite small sample size, depositional trends at the Haynie site generally align with those maintained across all sites in the Crow Canyon archaeofaunal database. Overall, dedicatory offerings are most abundant, followed by simple deposition and ceremonial discard. Unlike at Haynie, where NISP values for dedicatory offerings are significantly inflated and contribute less than 75% to the overall sample, no one type of deposition dominates across the entire Crow Canyon database (Figure 9).

Figure 9: Spatial distribution of canid remains across all sites included in the Crow Canyon archaeofaunal database. Relative abundance of canid specimens (%NISP) per depositional category are presented to visualize inter-site trends. Raw NISP counts are displayed beside each bar.

It is possible that some instances of canid deposition cannot be classified, or have been inappropriately classified according to the aforementioned three categories. Given the sheer volume and diversity of canine specimens in the Crow Canyon database, it is necessary to address the potential pitfalls of Hill’s interpretive approach. Some zooarchaeological research calls into question the breadth and scope of Hill’s descriptive frameworks. Monagle and Jones (2020) applied these categories to an analysis of Puebloan canid deposition, and found the archaeological data are contradicted by ethnographic accounts which blur the lines between ceremonial trash and dedicatory offerings.

Hill (2000) suggests that dedicatory offerings are differentiated from ceremonial trash based upon the function of the individual in death; dogs who have been ritually dispatched and deposited of as ceremonial trash serve no further role beyond that which they play in life, whereas dogs interred as dedicatory offerings continue to function as social and spiritual actors in death (Hill 2000). On the other hand, Pueblo groups consider discard a integral stage in the life of an object, welcoming its transition from one role in society to the next (Fladd et al. 2021; Walker and Berryman 2023). Hopi people, for instance, place particular emphasis on the intentional retirement of objects with important personal or ancestral connections (Walker and Berryman 2023).

The wolf mandible recovered in year one of archaeofaunal analysis (Dombrosky and Gilmore 2023) represents one such intersection between ceremonial discard and dedicatory offerings. The specimen is classified as simple deposition according to Figure 9, despite its recovery from surface room fill. The mandible cannot be classified as an element of ceremonial discard nor dedicatory offering according to Hill’s classification criteria (Hill 2000). The specimen nonetheless represents an intriguing instance of canid deposition. The wolf’s mandible may illustrate a specialized process of deposition by its very presence.

Pueblo ideas of canid deposition, both current and historic, offer an important avenue for unraveling the nature of human-canid interaction at Haynie. Moving forward, it is crucial to address the problematic nature of the term “trash” when discussing animal deposition in Pueblo society, as this terminology suggests a sense of indifference toward discarded materials following their deposition (Walker and Berryman 2023). The deposition of canid remains in midden deposits can reflect a range of personal and collective choices, enacted as a way to express, preserve, and commemorate key aspects of the living community (Fladd et al. 2021). We have therefore chosen to employ the term “ceremonial discard” in this analysis, and suggest that this term replace original terminology introduced by Hill (2000) to characterize discard in midden deposits.

Given these perspectives, it is apparent that dog deposition at the Haynie site is highly complex; colored by a vast and nuanced system of human-canid interaction. The application of frameworks such as Erica Hill’s (Hill 2000), which employ general descriptive categories, are useful for broader understandings of canine life in Pueblo society. For future analysis, it is crucial that Indigenous perspectives are prioritized as the primary interpretive framework for understanding canine life in the Mesa Verde region.

3.1.7 Large Mammal (n = 16)

The large mammal identification includes mammals larger than deer, and it includes specimens lacking morphological features required for more specific taxonomic levels. It is likely that it could incorporate large artiodactyls like elk and bison specimens. It should be incorporated in sensitivity analyses that rigorously assess conclusions based on these animals.

3.2 Aves (n = 219)

3.2.1 Large Birds (n = 135)

Birds larger than a mallard are considered large birds. This identification group is most likely dominated by Turkeys (Meleagris gallopavo) as they are one of the most frequent birds recovered from ancestral Pueblo sites. A high proportion of Turkey specimens are assigned to this group considering that there is considerable skeletal morphological overlap between Sandhill Crane (Grus canadensis) and Turkey (Hargrave and Emslie 1979). The Crow Canyon comparative collection does not, as of yet, include Sandhill Crane skeletal material. As a result, this group is frequently used as it represents a conservative identificaiton.

3.2.2 Galliformes (n = 56)

All Galliformes specimens identified so far in the Haynie assemblage are Turkey, but it is still surprising how few Turkey specimens there seem to be. For instance, Spielmann and Angstadt-Leto (1996) proposed the Turkey Index, which is the ratio of Turkeys to lagomorphs in an assemblage. The current Turkey Index for Haynie is 0.07, which is notably low. However, Driver (2002) proposed the Modified Turkey Index where the large bird identification group is included in the calculation. The current Modified Turkey Index value for Haynie is 0.20. This value is within the normal range for Pueblo I/Pueblo II sites in the central Mesa Verde region (Badenhorst and Driver 2009), but it is the lowest Modified Turkey Index for large assemblages in the Crow Canyon Research Database (Figure 10). Also worth noting, there appears to be an abundance of turkey specimens yet to be identified. It is doubtful that this low index value will remain constant as the faunal sample increases. How Turkey husbandry was managed at the community-level is an essential future area of research for the Northern Chaco Outliers Project. This line of inquiry is one way to delve deeper into aspects of cooperation and identity at the Lakeview Community.

Figure 10: The Modified Turkey Index across all faunal assemblages exceeding 1000 Number of Identified Specimens (NISP) in the Crow Canyon database. Currently, the Haynie site has the lowest Modified Turkey Index of all the central Mesa Verde assemblages (Tsama Pueblo is a Pueblo IV Northern Rio Grande site).

3.2.3 Medium Birds (n = 14)

Medium birds are considered larger than a robin and the size of a mallard or smaller. It is difficult to attribute the majority of this identification group to a single taxon, as it contains a variety of difficult-to-identify fragmented skeletal parts that could belong to numerous taxa.

3.2.4 Strigiformes (n = 5)

The owl specimens identified at the Haynie site are a terminal phalanx, second phalanx, first phalanx, femur, and ulna. Owl feathers are known to have been incorporated into dance paraphernalia and prayer sticks (Ladd 1963). It is, however, important to keep in mind that owls can be active taphonomic agents. Luckily, their signatures are well-known and include the presence of pellets, small mammal remains with little to no fragmentation, and visible signs of digestion on specimens (Andrews and Cook 1990; Fernández-Jalvo and Andrews 2016). Owls do not appear to be a taphonomic agent of concern at Haynie (Section 4).

3.2.5 Passeriformes (n = 3)

Three passeriform specimens have been identified so far in the Haynie assemblage: one small fragmented humerus generally identified to the order-level, a tibia fragment identified to Corvidae (the family comprising jays and crows), and one large carpometacarpus identified as Raven (Corvus corax).

3.2.6 Columbiformes (n = 2)

Pigeons and doves are in the order Columbiformes, and one ulna fragment and one sternum fragment were recovered from Haynie. These specimens are most likely a Mourning Dove (Zenaida macroura). An order-level identification was used because the Crow Canyon comparative collection does not include a Band-tailed Pigeon (Columba fasciata), which is the only other member in the order Columbiformes to consider in the region.

3.2.7 Accipitriformes (n = 3)

Two specimens were identified as general members of the order Accipitriformes: one general foot phalanx and one terminal phalanx that compared favorably to a Turkey Vulture (Cathartes aura). The other specimen was a distal tibiotarsus fragment identified to Golden Eagle (Aquila chrysaetos). This specimen was taken to the Museum of Southwestern Biology’s Division of Birds and identified with their skeletal comparative collection. It was distinguished from Bald Eagle (Haliaeetus leucocephalus) based off the morphology of the supratendinal bridge, which is more convex and more proximally robust in Golden Eagles compared to Bald Eagles. Eagles are acutely significant in Pueblo culture and their feathers are extremely valued (Beaglehole 1936; LaZar and Dombrosky 2022; Tyler 1991).

3.2.8 Anseriformes (n = 1)

One complete carpometacarpus was identified to the genus Anas spp. (mallards and other relatives) from the Haynie site. While this specimen could likely represent the common mallard (Anas platyrhynchos), Hargrave (1965) mentions the likely presence of Mexican Duck (Anas diazi) in the Four Corners region in the past. He based this assertion on the recovery of Mexican Duck feathers associated with Tusayan Polychrome sherds from Glen Canyon. The presence of different mallards and ducks is of potential biogeographic significance and the different members of this genus represent difficult-to-separate taxa.

3.3 Actinopterygii (n = 6)

This taxonomic class includes the ray-finned fishes. Six specimens have been identified to this general class: two ribs and three fragmented vertebrae centra (of which two of these specimens can be identified to the order-level, see below). Following Nelson (2006, 35), fish specimens should no longer be referred to as pisces, as it is an antiquated taxonomic term. Similarly, for fishes of inland North America, the use of osteichthyes should no longer be used (Nelson 2006, 83). This is so for two interrelated reasons. First, this term has been replaced by the Euteleostoma designation. It successfully describes a monophyletic clade that includes tetrapods. Secondly, since Euteloestoma includes tetrapods, it includes lobe-finned fishes (Sarcopterygii). Lobe-finned fishes—like the coelacanth (Actinistia)—are not native fishes in inland North America during the late Holocene (Cloutier and Forey 1991). The use of osteichthyes should be accordingly abandoned. Instead, Actinoptergygii (ray-finned fishes) should be used because it is a more accurate class-level designation for archaeofaunas from the U.S. Southwest/Mexican Northwest.

3.3.1 Cypriniformes (n = 2)

This order includes carps, minnows, and suckers, which are common fishes in the aridland streams of the U.S. Southwest (Minckley and Marsh 2009; Sublette et al. 1990). These two specimens are small intact vertebrae. The lateral ridge morphology of centra can be used to identify vertebrae of fishes from U.S. Southwestern archaeofaunas to the order-level. Common orders of fishes found in rivers in the U.S. Southwest include Cypriniformes, Siluriformes, Lepisosteiformes, and Salmoniformes, and each of these orders have distinct vertebral morphology. These specimens are also notably small. It is possible inhabitants of the Haynie site used non-targeted methods to capture fishes, such as seining (Dombrosky et al. 2022). A focus on fishing practices offers a basic way to understand aquatic habitat use associated with the Simon Draw watershed.

3.4 Reptilia (n = 1)

3.4.1 Squamata (n = 1)

One colubrid (family Colubridae) vertebra was identified from Haynie. This specimen was identified to the family-level using Olsen (1964). It likely represents a bullsnake (Pituophis catenifer sayi), which is one of the most common colubrids in Cortez, CO.

3.5 Oviparous Animal (n = 1)

Eggshell specimens were identified as oviparous animal. These specimens are likely from Turkey, but a general identification was given because a scanning electron microscope could not be used to assess mammillary cone morphology (Beacham and Durand 2007; Conrad et al. 2016; Lapham et al. 2016). It is possible, though unlikely, that this eggshell could be from a lizard.

4 Taphonomy

The most important taphonomic question surrounding the Haynie archaeofauna is: what is causing such a high degree of unidentifiability (Figure 1)? Last year, we posed the hypothesis that the high unidentifiability rate might be related to a high degree of artiodactyl exploitation (Dombrosky and Gilmore 2023). It could be that Haynie site occupants not only obtained a larger than average amount of artiodactyls (Figure 4), but that they intensively and extensively processed artiodactyl skeletal elements for within-bone nutrients as well (sensu Wolverton 2002). Haynie site occupants might have caused the high unidentifiability rate through their subsistence practices.

To support this hypothesis, we created a predictive model to test whether or not combinations of taphonomic variables could predict identifiable or unidentifiable specimens. Using logistic regression, we could confidently predict whether specimens were unidentifiable. The model we created indicated that one of the most important variables in predicting unidentifiability was whether a specimen had thick cortical bone. Considering that long bone shaft fragments have the highest incidence of thick cortical bone in our sample, we reasoned that these specimens are most likely linked to the most abundant medium-sized mammals recovered from the Haynie site (i.e., artiodactyls).

However, a crucial variable was left out of the predictive model we created: whether or not specimens come from areas of recent disturbance. This variable is crucial because the Haynie site exhibits extensive evidence for recently disturbed deposits, whether it be through mechanical disturbance, looter’s pits, or yard modification (Throgmorton et al. 2022, 2023). Unidentifiability could be linked to a high degree of artiodactyl fragmentation caused by recent disturbance. Here, we explore the relationship between unidentifiability through simple data visualization and hypothesis testing, we then add a new disturbance variable to the logistic regression model we developed last year. We also go a step further and compare our results to a different model type called random forest. This helps answer an additional question: could we develop more predictive accuracy—and clearer variable importance estimations—from using different model engines?

4.1 Identifiability by Recently Disturbed Contexts

There is a surprisingly similar ratio of unidentifiable to identifiable specimens from across disturbance contexts at Haynie given that almost half of the analyzed specimens are from contexts of recent disturbance (Figure 11). Approximately 78% of the specimens from recently disturbed contexts are unidentifiable, and about 67% of the specimens from undisturbed contexts are unidentifiable. Is this 11% difference significant given the sample sizes in each category? To test this question, we conducted a Pearson’s chi-square test of independence. It appears there is an association between identifiability and disturbance context (χ2 = 111.86, p < 0.01), with a low effect size (Cramér’s V = 0.12).

Inspecting the residuals (the difference between expected and observed counts) reveals an interesting pattern. There are larger residuals for the number of identifiable specimens in both disturbed and undisturbed contexts, while there are smaller residuals for unidentifiable specimens in each context (Figure 12). This suggests that unidentifiable specimens do not deviate from expected counts in each of the disturbance contexts.

Given the large sample size and the inspection of residuals, the association between identifiability and disturbance contexts is of very little practical significance. There is a non-random association, but it is of little magnitude or weight (see Wolverton et al. 2016). These results indicate that recently disturbed contexts will likely have little influence on predicting whether specimens are unidentifiable.

Figure 11: Percent of specimens that are unidentifiable and identifiable between recently disturbed and undisturbed contexts. Raw Number of Specimen (NSP) is reported within each bar.

Figure 12: Residuals (difference between expected and observed counts in each category) of identifiability between disturbance contexts. These residuals are an interesting reflection of one another. Residuals are similarly large (though in opposing directions) for identifiable specimens across contexts, and unidentifiable specimens have similarly moderate residuals (though again in opposing directions).

4.2 Updated Taphonomic Model

4.2.1 The Penalized Logistic Regression Model

Similar to last year, we randomly split our data into a training and test set, which comprises 80% and 20% of the data respectively. Our validation set consists of 80% of the data from the training set. This year, we stratified our training, test, and validation sets so that equal portions of identifiable and unidentifiable specimens are within each group. We compare penalized logistic regression and random forest models to predict unidentifiable specimens from this year’s sample. Dombrosky and Gilmore (2023) details the logic of creating a predictive model to understand taphonomic variable importance, we do not provide further justification here. The models are supplied with 17 predictor variables for every specimen: if it 1) came from an area of recent disturbance, if it had 2) thick cortical bone, 3) excavation damage, 4) carnivore damage, 5) at least one intact end, 6) a spiral fracture, 7) a transverse fracture, 8) an irregular break, if it was 9) a shaft fragment, 10) made into an artifact, 11) eroded, 12) gnawed by rodents, 13) splintered, 14) root etched, 15) burned, 16) who the analyst was, and 17) its maximum length.

The penalty term in the logistic regression model safeguards against highly correlated predictor variables, which can cause poor model performance. We used grid search on the validation set to tune 40 candidate penalty values, and we picked the one with the highest area under the Receiver Operator Curve (ROC) for our final model (Figure 13). The subsequent model was 87.37% accurate on the training set. Like last year, there was a large discrepancy between two other accuracy metrics. The Matthews correlation coefficient for the trained model is 0.68 while its F1 metric is 0.92. The F1 metric is markedly higher—indicating a model with high performance—because it describes how an event of interest is predicted. In this case, the event of interest is how well the model predicts unidentifiable specimens rather than identifiable ones (Figure 14). The model’s accuracy on our test set is remarkably similar, at 87.36%. As such, the logistic regression model fulfills the purpose of this analysis.

Figure 13: Area under the ROC (Receiver Operator Curve) for 40 penalized logistic regression models fit to our validation test set. The highest area under the ROC is highlighted with a blue dashed line.

Figure 14: Confusion matrix for the trained logistic regression model. This model accurately predicts unidentifiable specimens, but it performs poorly with identifiable specimens. This model was retained because we are most interested in what is driving the number of unidentifiable specimens.

4.2.2 Random Forest Comparison

Though the penalized regression model is generally high performing, it is possible that other predictive model types could provide higher accuracy and performance in different ways. We compare the logistic regression model to a random forest. Random forest models sometimes perform better than logistic regression models as the number of predictor variables increase in a dataset (Kirasich et al. 2018). Our model consists of a number of variables of unknown quality for prediction (i.e., noise versus explanatory variables), which underscores the need for optimal tuning parameters and model performance comparisons.

We tune two parameters in our random forest model: number of predictors randomly sampled at each tree and the minimum number of data points in a node required for further splitting. Parameters were selected using the highest area under the ROC (Figure 15). The tuned random forest model is 87.79% accurate on the training set, and there is a large discrepancy between the Mathews correlation coefficient (0.69) and F1 metric (0.92). These classification metrics for the trained random forest model are almost identical to those from logistic regression. And, indeed, the results are almost identical when comparing the two models’ performance on the test set (Figure 16). The random forest model does not increase or decrease predictive accuracy when compared to the logistic regression model.

Figure 15: Area under the ROC (Receiver Operator Curve) for two different tuning parameters (number of randomly selected predictors and minimal node size) for 40 different random forest models. Blue points represent the selected tuning values.

Figure 16: Receiver Operator Curve (ROC) comparing the tuned logistic regression and random forest models’ performance on the test set. Sensitivity is the true positive rate, which is how well the models classify truly unidentifiable specimens. The false positive rate (1 - specificity) is the probability that an identifiable specimen will get classified as an unidentifiable specimen. The dashed line represents the area at which a model will have random performance. Area Under the Curve (AUC) is a metric for model performance, where a model with an AUC of 1 equals perfect performance. Both the random forest and the logisitc regression models have an AUC of 88; they are nearly identical.

The random forest model does, however, provide a different view of variable importance, because it uses a series of drastically different steps to predict whether a specimen is unidentifiable or not (Figure 17). Disturbance context is the 7th most important variable in the logistic regression model. The two most important variables, which detract from unidentifiability, are whether or not a specimen has one intact end and the specimen’s length. The two most important variables adding to unidentifiability are whether the specimen has thick cortical bone or whether the specimen has a spiral fracture. For the random forest model, disturbance context is the 10th most important variable. The top three most important variables are whether or not there is an intact end, the length of the specimen, and/or whether it is a shaft fragment.

Figure 17: Variable importance plot for the unidentifiability model developed here (Greenwell and Boehmke 2020). The sign of the coefficient is plotted next to each bar for the logistic regression model, which indicates whether the variable adds or subtracts from unidentifiability (the event of interest). The blue bar highlights where recently distrubed contexts fall out in predicting unidentifiability. Disturbance does not significantly impact unidentifiability.

The most important variable across both models is whether or not a specimen has at least one intact end. The more intact a specimen is the more identifiable it is. We believe this result is a sign of high data quality. An assemblage with high identifiability on specimens lacking morphologically distinct features (i.e., intact ends) would be alarming.

Our new models further support the hypothesis that medium-to-large mammal fragmentation is a key taphonomic feature of this assemblage. However, the above results also indicate that disturbance context does not heavily influence unidentifiability. We believe artiodactyl exploitation mixed with intensive and extensive within bone nutrient processing could be the cause of the high unidentifiability rate.

High numbers of artiodactyl specimens relative to other staple hunted animals (lagomorphs) paired with high fragmentation pose a potential conundrum. On one hand, these specimens could indicate efficient foraging practices, where Haynie hunters were able to consistently acquire high-ranked prey. On the other hand, foraging efficiency could have been low since large game were so thoroughly processed. The calculation of fragmentation rates per skeletal part and across the skeleton of individual artiodactyls will help resolve this issue (sensu Wolverton 2002). Understanding these patterns through time will be especially important to consider. Future effort should be focused on assessing whether artiodactyl fragmentation increases through time or if the amount of fragmentation remains stable through Haynie’s occupational sequence. It could be that culinary practices were thorough—requiring high levels of processing—and remained stable through time.

This future analysis will also help address a glaring and basic issue: is fragmentation artificially inflating NISP values, obscuring the relatively high artiodactyl index at Haynie? Next year’s report will focus explicitly on this issue. We will calculate the Minimum Number of Individuals (MNI), Minimum Number of Elements (MNE), %whole, and NISP:MNE at Haynie and across sites in the Crow Canyon Research Database.

5 Taxonomic Diversity and Representative Sampling

Last year, we wanted to understand whether or not Haynie had a more diverse fauna compared to other sites. The Lakeview Community is located in a natural corridor (Simon Draw) that could have facilitated the movement of both birds and large game, supplying highly diverse wild resources to Haynie site residents. We compared rarefaction curves across three other archaeological sites to test this hypothesis (Shields Pueblo, Sand Canyon Pueblo, and Ponsipa’akeri). Rarefaction curves allow us to compare taxonomic diversity across multiple sites (while controlling for sample size) and to assess whether sampling efforts have been sufficient enough to provide accurate taxonomic representation. In our initial analysis we found that Haynie was not remarkably diverse and that sampling still was not suffecient enough to provide an accurate estimation of taxonomic diversity. Here, we update the rarefaction curves with sites across the entire Crow Canyon Research Database that exceed 2000 NISP.

This analysis relies on a non-standard quantitative unit to estimate taxonomic richness called Unique Identifications Types (UITs). It is a tally of the different identification types present in a specific context, meaning it can include standard taxonomic identifications (e.g., Odocoileus sp.) and non-standard identifications (e.g., medium artiodactyl). This unit serves as a proxy for taxonomic richness to help gauge patterns in sampling and recovery. We prefer this unit over the common Number of Taxa (NTAXA) for three reasons: 1) it is simpler to calculate when dealing with large mixed assemblages identified to a variety of taxonomic levels, 2) it does not require the selection of an arbitrary taxonomic group from which to aggregate all lower units within, and 3) it is strongly correlated with NTAXA when calculated at multiple levels of taxonomic resolution (Dombrosky and Gilmore 2023, fig. 18).

Taxonomic diversity accumulates fastest within the Haynie site fauna compared to six other large assemblages from the central Mesa Verde region (Figure 18). However, we do not yet consider the rate of accumulation significantly different enough to argue that species diversity at Haynie is somehow anomalous. One new identification type is added every 116 NISP at Haynie, while the other sites usually add a new identification type about every 162 NISP. Future work will focus on factoring in the duration of site occupation at each site (sensu Varien and Potter 1997; Varien and Mills 1997; Varien and Ortman 2005). Controlling for site occupation and increasing sample size will clarify these patterns. Additionally, we view analyzing 116 identifiable specimens for every new identification type low sampling effort. Based on these findings, we also conclude that sampling effort has not yet been sufficient enough to accurately capture all of the taxonomic diversity at Haynie. The final word on taxonomic diversity will have to wait until we analyze more specimens. If the rate of change continues its upward trajectory, and does not level-off, Haynie may indeed be interpreted as anomalously diverse faunal assemblage in the central Mesa Verde region.

Figure 18: Rarefied species accumulation curves. We randomly selected 2000 identified specimens from each archaeofaunal assemblage and replicated this 500 times, then calculated the mean Unique Identification Types (UITs) as specimens accumulate, along with standard deviation (the gray ribbon). We calculated the Average Rate of Change (ARC) for each curve from 1500 to 2000 NISP, which serves to inidicate whether sampling efforts have been sufficient at this point in the analysis. ARC is represented by the slope of the thick, blue dashed line.

6 Conclusion

Archaeofaunal analysis at the Haynie site continues to point to new and exciting areas of research. Inter- and intrasite analysis of canid deposition reveals a complex system of human-canine interaction at Haynie and across the central Mesa Verde region, challenging conventional conceptualizations of domestication. Further analysis will consider canine age profiles and associations with non-canid taxa, painting a more vivid picture of canine-human relationality in Ancestral Pueblo society. Haynie also currently has one of the highest artiodactyl index values in the region, and some of the lowest lagomorph and turkey indices. Taphonomic analysis continues to suggest that high numbers of artiodactyls and the processing of their skeletal parts is a primary cause of the high unidentifiability rate at the site, even when recently disturbed contexts are considered. Next year’s report will explicitly include an in-depth analysis of fragmentation patterns at Haynie and across assemblages in the Crow Canyon Research Database. Taxonomic diversity showed that the Haynie archaeofaunal assemblage is not significantly diverse compared to other sites, but that it is likely the most diverse faunal assemblage in the Crow Canyon database. If the rate of accumulation for new identification types continues, then it is possible the Haynie site could be one of the most diverse faunal assemblages in the central Mesa Verde region. Continued sampling will help clarify this possibility. Zooarchaeology is a crucial component of the Northern Chaco Outliers Project and will continue to provide exciting information to the archaeology of human-environment interaction in the U.S. Southwest and beyond.

7 Acknowledgements

Thanks to Mac Marston, Susan Ryan, Dave Satterwhite, and Catherine West for helpful comments and edits on previous drafts of this report. Many continued thanks to Steve Wolverton for his zooarchaeological support and inspiration.

8 References

Adams, Bradley J., and Pam J. Crabtree. 2012. Comparative Osteology a Laboratory and Field Guide of Common North American Animals. Boston, MA: Elsevier/Academic Press.
Ambrose, Stanley H. 1990. “Preparation and Characterization of Bone and Tooth Collagen for Isotopic Analysis.” Journal of Archaeological Science 17 (4): 431–51. https://doi.org/10.1016/0305-4403(90)90007-R.
Andrews, Peter, and Jill Cook. 1990. Owls, Caves and Fossils: Predation, Preservation and Accumulation of Small Mammal Bones in Caves, with an Analysis of the Pleistocene Cave Faunas from Westbury-Sub-Mendip, Somerset, UK. Chicago, IL: University of Chicago Press.
Badenhorst, Shaw, and Jonathan C. Driver. 2009. Faunal Changes in Farming Communities from Basketmaker II to Pueblo III (A.D. 1-1300) in the San Juan Basin of the American Southwest.” Journal of Archaeological Science 36 (9): 1832–41. https://doi.org/10.1016/j.jas.2009.04.006.
Badenhorst, Shaw, Jonathan C. Driver, and Steve Wolverton. 2023. The Exploitation of Rodents in the Mesa Verde Region.” In Research, Education, and American Indian Parternships at the Crow Canyon Archaeological Center, edited by Susan C. Ryan, 325–34. Salt Lake City: The University of Utah Press. https://doi.org/10.5876/9781646424597.c020.
Baker, Brenda J., Dupras Tosha L., and Matthew W. Tocheri. 2005. The Osteology of Infants and Children. Vol. 12. College Station, TX: Texas A&M University Press.
Baker, Bruce W., and Edward P. Hill. 2003. Beaver (Castor Canadensis).” In Wild Mammals of North America: Biology, Management, and Conservation, edited by George A. Feldhamer, Bruce C. Thompson, and Joseph A. Chapman, 2nd ed., 288–310. Baltimore, MD: The Johns Hopkins University Press.
Beacham, E. Bradley, and Stephen R. Durand. 2007. Eggshell and the Archaeological Record: New Insights into Turkey Husbandry in the American Southwest.” Journal of Archaeological Science 34 (10): 1610–21. https://doi.org/10.1016/j.jas.2006.11.015.
Beaglehole, Ernest. 1936. Hopi Hunting and Ritual. Yale University Publications in Anthropology 4. New Haven, CT: Yale University Press.
Bocinsky, R. Kyle, Jason A. Cowan, Timothy A. Kohler, and C. David Johnson. 2012. How Hunting Changes the VEP World, and How the VEP World Changes Hunting.” In Emergence and Collapse of Early Villages: Models of Central Mesa Verde Archaeology, edited by Timothy A. Kohler and Mark Varien, 145–52. Berkeley: University of California Press.
Broughton, Jack M., Michael D. Cannon, Frank E. Bayham, and David A. Byers. 2011. Prey Body Size and Ranking in Zooarchaeology: Theory, Empirical Evidence, and Applications from the Northern Great Basin.” American Antiquity 76 (3): 403–28. https://doi.org/10.7183/0002-7316.76.3.403.
Broughton, Jack M., and Shawn D. Miller. 2016. Zooarchaeology and Field Ecology: A Photographic Atlas. Salt Lake City: The University of Utah Press.
Center, Crow Canyon Archaeological. 2001. Crow Canyon Archaeological Center Field Manual. https://crowcanyon.org/ResearchReports/FieldManual/hold_Field_Manual_all.pdf.
Chavez, Angela Marie. 2008. Comparative Vertebral Morphology in Medium-Sized North American Artiodactyla.” The Artifact 46: 1–17.
Cloutier, Richard, and Peter L. Forey. 1991. Diversity of Extinct and Living Actinistian Fishes (Sarcopterygii).” In The Biology of Latimeria Chalumnae and Evolution of Coelacanths, edited by John A. Musick, Michael N. Bruton, and Eugene K. Balon, 59–74. Springer.
Cohen, Alan, and D. Serjeantson. 1996. A Manual for the Identification of Bird Bones from Archaeological Sites. London: Archetype Publications.
Cole, Kasey E., Peter M. Yaworsky, and Isaac A. Hart. 2022. “Evaluating Statistical Models for Establishing Morphometric Taxonomic Identifications and a New Approach Using Random Forest.” Journal of Archaeological Science 143: 105610. https://doi.org/10.1016/j.jas.2022.105610.
Conrad, Cyler, Emily Lena Jones, Seth D. Newsome, and Douglas W. Schwartz. 2016. One Isotopes, Eggshell and Turkey Husbandry at Arroyo Hondo Pueblo.” Journal of Archaeological Science: Reports 10: 566–74. https://doi.org/10.1016/j.jasrep.2016.06.016.
Dombrosky, Jonathan, Caitlin S. Ainsworth, Abigail A. Judkins, Jana Valesca Meyer, Michael A. Adler, and Emily Lena Jones. 2023. “Reconsidering Archaeological Garden Hunting: A View from the Northern U.S. Southwest.” Journal of Archaeological Science: Reports 51: 104194. https://doi.org/10.1016/j.jasrep.2023.104194.
Dombrosky, Jonathan, and Eric Gilmore. 2023. Year One of Archaeofaunal Analysis at the Haynie Site (5MT1905).” In Excavation and Additional Studies at The Haynie Site (5MT1905) by the Crow Canyon Archaeological Center Annual Report 2022, edited by Kellam Throgmorton, Susan C. Ryan, Jonathan Dombrosky, Benjamin Bellorado, Steven Copeland, Kate Hughes, and Jamie Merewether, 45–88. https://www.crowcanyon.org/wp-content/uploads/2023/02/ncop-annual-report-2022-final.pdf.
Dombrosky, Jonathan, Thomas F. Turner, Alexandra Harris, and Emily Lena Jones. 2022. Body Size from Unconventional Specimens: A 3D Geometric Morphometrics Approach to Fishes from Ancestral Pueblo Contexts.” Journal of Archaeological Science 142: 105600. https://doi.org/10.1016/j.jas.2022.105600.
Driver, Jonathan C. 1992. Identification, Classification and Zooarchaeology.” Circaea 9 (1): 35–47.
———. 2002. Faunal Variation and Change in the Northern San Juan Region.” In Seeking the Center Place: Archaeology and Ancient Communities in the Mesa Verde Region, edited by Mark D. Varien and Richard H. Wilshusen, 143–60. Salt Lake City: The Univeristy of Utah Press.
———. 2006. Crow Canyon Archaeological Center Manual for Description of Vertebrate Remains. 8th ed. Cortez, CO: Crow Canyon Archaeological Center.
———. 2011. Human Impacts on Animal Populations in the American Southwest.” In Movement, Connectivity and Landscape Change in the Ancient Southwest, edited by Margaret C. Nelson and Colleen Strawhacker, 179–98. Boulder: University Press of Colorado.
———. 2011. Identification, Classification and Zooarchaeology (Featured Reprint and Invited Comments).” Ethnobiology Letters 2: 19–39. https://doi.org/10.14237/ebl.2.2011.32.
Driver, Jonathan C., and Joshua R. Woiderski. 2008. Interpretation of the "Lagomorph Index" in the American Southwest.” Quaternary International 185: 3–11. https://doi.org/10.1016/j.quaint.2007.09.022.
Elbroch, Mark. 2006. Animal Skulls: A Guide to North American Species. Mechanicsburg, PA: Stackpole Books.
Ellyson, Laura. 2014. Resource Intensification of Small Game Use at Goodman Point, Southwestern Colorado.” Master's Thesis, Denton: Department of Geography, University of North Texas. https://digital.library.unt.edu/ark:/67531/metadc699883/.
Fernández-Jalvo, Yolanda, and Peter Andrews. 2016. Atlas of Taphonomic Identifications: 1001+ Images of Fossil and Recent Mammal Bone Modification. Springer.
Fladd, Samantha G, Saul L Hedquist, and E Charles Adams. 2021. “Trash Reconsidered: A Relational Approach to Deposition in the Pueblo Southwest.” Journal of Anthropological Archaeology 61: 101268.
France, Diane L. 2009. Human and Nonhuman Bone Identification: A Color Atlas. Boca Raton, FL: CRC Press.
Gilbert, B. Miles. 1980. Mammalian Osteology. Laramie, WY: B. M. Gilbert.
Gilbert, B. Miles, Larry D. Martin, and Howard G. Savage. 1981. Avian Osteology. Laramie, WY: B. M. Gilbert.
Gobalet, Kenneth W. 2001. A Critique of Faunal Analysis: Inconsistency Among Experts in Blind Tests.” Journal of Archaeological Science 28 (4): 377–86. https://doi.org/10.1006/jasc.2000.0564.
Grayson, Donald K. 1984. Quantitative Zooarchaeology: Topics in the Analysis of Archaeological Faunas. Orlando, FL: Academic Press. https://www.sciencedirect.com/book/9780122972805/quantitative-zooarchaeology.
Grayson, Donald K, and Carol J Frey. 2004. “Measuring Skeletal Part Representation in Archaeological Faunas.” Journal of Taphonomy 2 (1): 27–42.
Greenwell, Brandon M., and Bradley C. Boehmke. 2020. Variable Importance Plots—An Introduction to the vip Package.” The R Journal 12 (1): 343–66. https://doi.org/10.32614/RJ-2020-013.
Guiry, Eric, Trevor J. Orchard, Suzanne Needs-Howarth, and Paul Szpak. 2021. Isotopic Evidence for Garden Hunting and Resource Depression in the Late Woodland of Northeastern North America.” American Antiquity 86 (1): 90–110. https://doi.org/10.1017/aaq.2020.86.
Hargrave, Lyndon L. 1965. Archaeological Bird Bones from Chapin Mesa, Mesa Verde National Park.” American Antiquity 31 (2): 156–60. https://doi.org/10.1017/S0081130000004494.
Hargrave, Lyndon L., and Steven D. Emslie. 1979. Osteological Identification of Sandhill Crane Versus Turkey.” American Antiquity 44 (2): 295–99. https://doi.org/10.2307/279079.
Hastie, Trevor, Robert Tibshirani, and Jerome H. Friedman. 2009. The Elements of Statistical Learning: Data Mining, Inference, and Prediction. Springer.
Hill, Erica. 2000. The Contextual Analysis of Animal Interments and Ritual Practice in Southwestern North America.” KIVA 65 (4): 361–98. https://doi.org/10.1080/00231940.2000.11758417.
Hillson, Simon. 1986. Teeth. New York: Cambridge University Press.
———. 1996. Mammal Bones and Teeth: An Introductory Guide to Methods of Identification. London: Institute of Archaeology, University College London.
Jacobson, Jodi A. 2003. Identification of Mule Deer (Odocoileus Hemionus) and White-Tailed Deer (Odocoileus Virginianus) Postcranial Remains as a Means of Determining Human Subsistence Strategies.” Plains Anthropologist 48 (187): 287–97. https://doi.org/10.1080/2052546.2003.11949269.
James, Gareth, Daniela Witten, Trevor Hastie, and Robert Tibshirani. 2013. An Introduction to Statistical Learning: With Applications in R. Springer.
Kirasich, Kaitlin, Trace Smith, and Bivin Sadler. 2018. “Random Forest Vs Logistic Regression: Binary Classification for Heterogeneous Datasets.” SMU Data Science Review 1 (3). https://scholar.smu.edu/datasciencereview/vol1/iss3/9/.
Kuhn, Max, and Kjell Johnson. 2013. Applied Predictive Modeling. Springer.
Kuhn, Max, and Julia Silge. 2022. Tidy Modeling with R. Sebastopol, CA: O’Reilly Media, Inc. https://www.tmwr.org/.
Ladd, Edmund J. 1963. Zuni Ethno-Ornithology.” Master's Thesis, Albuquerque: Department of Anthropology, University of New Mexico.
Lapham, Heather A., Gary M. Feinman, and Linda M. Nicholas. 2016. Turkey Husbandry and Use in Oaxaca, Mexico: A Contextual Study of Turkey Remains and SEM Analysis of Eggshell from the Mitla Fortress.” Journal of Archaeological Science: Reports 10: 534–46. https://doi.org/10.1016/j.jasrep.2016.05.058.
LaZar, Miranda, and Jonathan Dombrosky. 2022. Exploring Raptor Exchange.” Archaeology Southwest Magazine 35 (1 & 2): 29–31. https://www.archaeologysouthwest.org/product/asw35-1-2/.
Lyman, R. Lee. 2002. Taxonomic Identification of Zooarchaeological Remains.” The Review of Archaeology 23 (2): 13–20.
———. 2008. Quantitative Paleozoology. Cambridge: Cambridge University Press.
———. 2019. Assumptions and Protocol of the Taxonomic Identification of Faunal Remains in Zooarchaeology: A North American Perspective.” Journal of Archaeological Method and Theory 26 (4): 1376–1438. https://doi.org/10.1007/s10816-019-09414-0.
Minckley, Wendell L., and Paul C. Marsh. 2009. Inland Fishes of the Greater Southwest: Chronicle of a Vanishing Biota. Tucson: The University of Arizona Press.
Monagle, Victoria, Cyler Conrad, and Emily Lena Jones. 2018. “What Makes a Dog? Stable Isotope Analysis and Human-Canid Relationships at Arroyo Hondo Pueblo.” Open Quaternary 4 (1).
Monagle, Victoria, and Emily Lena Jones. 2020. Dog Life and Death in an Ancestral Pueblo Landscape.” In Dogs: Archaeology beyond Domestication, edited by Brandi Bethke and Amanda Burtt, 45–71. Gainesville: University Press of Florida.
Nelson, Joseph S. 2006. Fishes of the World. Hoboken, NJ: John Wiley.
Olsen, Stanley John. 1964. Fish, Amphibian and Reptile Remains from Archaeological Sites. Cambridge, MA: Peabody Museum.
———. 1964. Mammal Remains from Archaeological Sites. Cambridge, MA: Peabody Museum.
———. 1979. Osteology for the Archaeologist. Cambridge, MA: Peabody Museum.
R Core Team. 2023. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing. https://www.R-project.org/.
Rawlings, Tiffany A. 2006. Faunal Analysis and Meat Procurement: Reconstructing the Sexual Division of Labor at Shields Pueblo, Colorado.” Doctoral Dissertation, Burnaby, BC: Department of Archaeology, Simon Fraser University. https://summit.sfu.ca/item/2374.
Roze, Uldis, and Linda M. Ilse. 2003. Porcupine (Erethizon Dorsatum).” In Wild Mammals of North America: Biology, Management, and Conservation, edited by George A. Feldhamer, Bruce C. Thompson, and Joseph A. Chapman, 2nd ed., 371–80. Baltimore, MD: The Johns Hopkins University Press.
Ryan, Susan C. 2016. The Northern Chaco Outliers Project: A Proposal to Conduct Archaeological Testing at the Haynie Site, Southwestern Colorado.” Manuscript on file Crow Canyon Archaeological Center, Cortez, CO. Submitted to Colorado State Historic Preservation Office, Denver, CO.
Semanko, Amanda, and Frank C Ramos. 2022. “Teaching an Old Dog New Tricks: Implications for Isotopic Studies of Southwest Dogs.” Kiva 88 (3): 327–46.
Simon, Rebecca L., Susan C. Ryan, Shanna R. Diederichs, Kari L. Schleher, Caitlin A. Sommer, Steven R. Copeland, and Grant D. Coffey. 2017. The Northern Chaco Outliers Project Annual Report, 2017 Field Season.” Cortez, CO: Crow Canyon Archaeological Center. https://crowcanyon.org/ResearchReports/NorthernChacoOutliers/interim_reports/NCOP_AnnualReport2017_FINAL.pdf.
Smart, Tamela S. 2009. Carpals and Tarsals of Mule Deer, Black Bear and Human: An Osteology Guide for the Archaeologist.” Master's Thesis, Pullman: Department of Anthropology, Western Washington University. https://cedar.wwu.edu/cgi/viewcontent.cgi?referer=&httpsredir=1&article=1018&context=wwuet.
Spielmann, Katherine A., and Eric A. Angstadt-Leto. 1996. Hunting, Gathering, and Health in the Prehistoric Southwest.” In Evolving Complexity and Environmental Risk in the Prehistoric Southwest, edited by Joseph A. Tainter, and Bonnie B. Tainter, 79–106. Reading, MA: CRC Press.
Sublette, James E., Michael D. Hatch, and Mary F. Sublette. 1990. The Fishes of New Mexico. Albuquerque: University of New Mexico Press.
Throgmorton, Kellam, Susan C. Ryan, Benjamin Bellorado, Jeremy Grundvig, Steve Copeland, and Timothy Wilcox. 2022. Excavation and Additional Studies at The Haynie Site (5MT1905) by the Crow Canyon Archaeological Center.” Annual Report 2021. Cortez, CO: Crow Canyon Archaeological Center. https://www.crowcanyon.org/wp-content/uploads/2022/02/ncop-annual-report-2021-final-1.pdf.
Throgmorton, Kellam, Susan C. Ryan, Jonathan Dombrosky, Benjamin Bellorado, Steven Copeland, Kate Hughes, and Jamie Merewether. 2023. Excavation and Additional Studies at The Haynie Site (5MT1905) by the Crow Canyon Archaeological Center.” Annual Report 2022. Cortez, CO: Crow Canyon Archaeological Center. https://crowcanyon.org/wp-content/uploads/2023/02/ncop-annual-report-2022-final.pdf.
Tyler, Hamilton A. 1991. Pueblo Birds and Myths. Flagstaff, AZ: Northland Publishing Company.
Varien, Mark D., and Barbara J. Mills. 1997. Accumulations Research: Problems and Prospects for Estimating Site Occupation Span.” Journal of Archaeological Method and Theory 4 (2): 141–91. https://doi.org/10.1007/BF02428057.
Varien, Mark D., and Scott G. Ortman. 2005. Accumulations Research in the Southwest United States: Middle-Range Theory for Big-Picture Problems.” World Archaeology 37 (1): 132–55. https://doi.org/10.1080/0043824042000329603.
Varien, Mark D., and James M. Potter. 1997. Unpacking the Discard Equation: Simulating the Accumulation of Artifacts in the Archaeological Record.” American Antiquity 62 (2): 194–213. https://doi.org/10.2307/282506.
Walker, William H, and Judy Berryman. 2023. “Ritual Closure: Rites de Passage and Apotropaic Magic in an Animate World.” Journal of Archaeological Method and Theory 30 (2): 449–94.
Welker, Martin H., David A. Byers, and Sarah B. McClure. 2021. ‘I Wanna Be Your Dog’: Evaluating the Efficacy of Univariate and Multivariate Methods for Differentiating Domestic and Wild Canids in North America.” International Journal of Osteoarchaeology 31 (2): 196–206. https://doi.org/10.1002/oa.2939.
White, Tim D., Michael T. Black, and Pieter A. Folkens. 2012. Human Osteology. 3rd ed. Cambridge, MA: Academic Press.
Wickham, Hadley. 2010. A Layered Grammar of Graphics.” Journal of Computational and Graphical Statistics 19 (1): 3–28. https://doi.org/10.1198/jcgs.2009.07098.
Wickham, Hadley, Mara Averick, Jennifer Bryan, Winston Chang, Lucy D’Agostino McGowan, Romain François, Garrett Grolemund, Alex Hayes, Lionel Henry, and Jim Hester. 2019. Welcome to the Tidyverse.” Journal of Open Source Software 4 (43): 1686. https://doi.org/10.21105/joss.01686.
Wolverton, Steve. 2002. NISP:MNE and %Whole in Analysis of Prehistoric Carcass Exploitation.” North American Archaeologist 23 (2): 85–100. https://doi.org/10.2190/egdq-cq1q-lld2-h3tp.
———. 2013. Data Quality in Zooarchaeological Faunal Identification.” Journal of Archaeological Method and Theory 20 (3): 381–96. https://doi.org/10.1007/s10816-012-9161-4.
Wolverton, Steve, Jonathan Dombrosky, and R. Lee Lyman. 2016. “Practical Significance: Ordinal Scale Data and Effect Size in Zooarchaeology.” International Journal of Osteoarchaeology 26 (2): 255–65. https://doi.org/10.1002/oa.2416.
Wolverton, Steve, and Lisa Nagaoka. 2018. Zooarcheology: Investigating Past Interactions Between Humans and Other Animals.” In Ethnozoology, edited by Rômulo Romeu Nóbrega Alves and Ulysses Paulino Albuquerque, 25–43. Cambridge, MA: Academic Press.

Citation

BibTeX citation:
@incollection{dombrosky2024,
  author = {Dombrosky, Jonathan and Feldstein, Ahna},
  editor = {David Satterwhite, R. and C. Ryan, Susan and Throgmorton,
    Kellam and Copeland, Steve and Hughes, Kate and Merewether, Jamie
    and Sinensky, Reuven},
  title = {Year {Two} of {Archaeofaunal} {Analysis} at the {Haynie}
    {Site} {(5MT1905)}},
  booktitle = {Excavation and Additional Studies at The Haynie Site
    (5MT1905) by Crow Canyon Archaeological Center Annual Progress
    Report 2023},
  pages = {70-104},
  date = {2024},
  url = {https://crowcanyon.org/wp-content/uploads/2024/02/NCOP-Annual-Report-2023_-Final.pdf},
  langid = {en}
}
For attribution, please cite this work as:
Dombrosky, Jonathan, and Ahna Feldstein. 2024. “Year Two of Archaeofaunal Analysis at the Haynie Site (5MT1905).” In Excavation and Additional Studies at The Haynie Site (5MT1905) by Crow Canyon Archaeological Center Annual Progress Report 2023, edited by R. David Satterwhite, Susan C. Ryan, Kellam Throgmorton, Steve Copeland, Kate Hughes, Jamie Merewether, and Reuven Sinensky, 70–104. https://crowcanyon.org/wp-content/uploads/2024/02/NCOP-Annual-Report-2023_-Final.pdf.